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Abstract 

The standard calculations of the synchrotron emission from charged particles in 
magnetic fields does not apply when the energy losses of the particles are so se- 
vere that their energy is appreciably degraded during one Larmor rotation. In these 
conditions, the intensity and spectrum of the emitted radiation depend on the ob- 
servation time Tobs- the standard result is recovered only in the limit Tq^s <C Tioss, 
where Tioss is the time for synchrotron losses. In this case the effects of the radi- 
ation backreaction cannot be detected by the observer. We calculate the emitted 
power of the radiation in the most general case, naturally including both the cases 
in which the backreaction is relevant and the standard case, where the usual result 
is recovered. Finally we propose several scenarios of astrophysical interest in which 
the effects of the backreaction cannot and should not be ignored. 



1 Introduction 



As pointed out in [1] (hereafter paper I), some aspects of the synchrotron 
emission of high energy particles did not receive proper attention in the lit- 
erature, mainly because the standard treatment proved to be valid over most 
of the energy range of interest for astrophysical applications. Recently, this 
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range of interest changed considerably, mainly due to the discovery of radia- 
tion processes at ultra-high energies. Nevertheless the important corrections 
to synchrotron emission in this regime have been ignored and the standard 
calculations have been adopted. 

In paper I we explored the generalization of the calculations of the synchrotron 
emission from an ensemble of particles radiating coherently. We considered 
there the cases of bunches of particles both in the monoenergetic case and in 
the case of a spectrum of particles, and for each we established the criteria for 
the synchrotron emission to be coherent. 

In the present paper we explore an effect that becomes relevant at sufficiently 
high energies, when the energy lost by a particle during one Larmor gyra- 
tion becomes comparable with the energy of the particle itself. We call this 
effect backreaction, with may be a slightly improper term. We find that in 
this regime there are important corrections to the spectra of the emitted radi- 
ation in numerous scenarios currently discussed in the literature, confirming 
but extending previous findings of Ref. [2], where the case of a monoener- 
getic distribution of radiating electrons was considered. In this paper we also 
study the situation of an ensemble of particles with arbitrary energy spectra 
and find an approximate analytical solution of the problem of the synchrotron 
backreaction, that may be useful to estimate the magnitude of the effect. The 
formalism used here, introduced in paper I and briefly summarized in this 
paper, allows us to take naturally into account possible coherence effects in 
the synchrotron radiation, together with the backreaction. Cases of astrophys- 
ical interest in which the effects of the backreaction are supposed to play an 
important role will be discussed. 

The paper is structured as follows: in section 2 we describe the backreaction 
and the effects that can be expected on simple basis. In section 3 we calculate 
the spectrum of the radiation emitted by monoenergetic particles. In section 
4 we generalize the calculation to the case of a spectrum of radiating parti- 
cles and we present analytical approximations that allow us to estimate the 
effects of the backreaction without (or before) being involved in the detailed 
calculations. We present our conclusions in section 5. 



2 The backreaction 

The rate of energy losses of a particle with mass m, Lorentz factor 7 and 
charge q in a, magnetic field B can be written in the well known form 



2q^ 
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where for simplicity we assumed that the electron moves perpendicular to the 
direction of the magnetic field. The time for a Larmor rotation can be easily 
calculated as tb = 2'nE/qBc. The (often) hidden assumption of the standard 
calculations of synchrotron radiation is that the energy lost by a particle in 
the time tb is negligible compared with the particle energy E. This condition 
is fulfilled when 

Tb ^ fSm^c^Y^^ 

— < 1 ^ < — — - 2 

where tr is the typical time of energy losses of the particle, defined by the 
expression 

J 1 dE 2q^B^ ^ 

^ E dt ^ 2,m'^c^ ^ ' 



For energies larger than the value found in eq. (2), the standard calculations 
fail and a new approach is needed. In the present paper we introduce this novel 
approach and show that the difference between the predicted spectra and the 
standard ones may be extremely important and in general cannot be ignored. 

Note that from equation (2) the region where the backreaction becomes im- 
portant is numerically given by 

E » 1.9 X 10^5^^.. GeV (4) 

for electrons, and by 

£; » 6.5 X 10^°Sg:/L GeV (5) 

for protons. Here the magnetic field Bcauss is in Gauss. A visual picture of 
the region in the plane E — B where the backreaction is relevant is provided 
in fig. 1 for both electrons and protons: the backreaction should be taken into 
account above the solid lines reported in fig. 1. 

3 Spectrum of monoenergetic particles 



In this section we compute the spectrum radiated by Z particles, all having the 
same initial Lorentz factor 70, as detected by an observer during an observation 
time Tohs (this calculation was also carried out with some differences in the 
approach, in [2]). This situation, as shown in [2] is equivalent to calculate 
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Fig. 1. Region in plane E — B where the backreaction becomes relevant. The case of 
electrons and protons is illustrated. 

the power radiated by a single particle with Lorentz factor 70, as detected 
by Z observers during the observation time. Although trivial for the standard 
scenario, this is not trivial for the situation of interest here. 

We also demonstrate that for observation times smaller than the time for 
appreciable losses, the standard result is recovered. 

A nice didactical approach illustrated in [2] allows to show that the result 
of the backreaction, that should be derived by the solution of the Lorentz- 
Dirac equation, in the end can be determined in the simple way of taking into 
account the time dependence of the Lorentz factor of the radiating particle. 

Based on eq. (1), we can write the time dependence of the Lorentz factor of a 
particle subject to synchrotron energy losses in the following form 



where ojr — ojkIo = T/j conveniently defines the time scale for the energy 
losses of a particle. In [2] it was shown that this expression keeps its validity 
also in case of strong backreaction, at least as long as quantum effects remain 
negligible. These effects enter the calculation only when the typical energy of 
the radiated photon becomes comparable to the energy of the particle, so that 
quantum recoil needs to be accounted for. We will not be involved here in 
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these extreme cases, although we behave that some interesting physics coTild 
be learned there. To our knowledge, no paper exists in the literature treating 
the problem of the synchrotron emission with both the backreaction and the 
quantum recoil at the same time. 

Based on the formalism illustrated in detail in paper I, we associate to each 
particle a phase factor a. Using the usual relativistic beaming condition, im- 
plying that significant synchrotron emission is confined in a narrow cone of 
aperture ~ I/7, we deduce that a given particle illuminates the observer at 
those retarded times that are related to the phase a by the equation 

tujB{t) - a (7) 



where LOsif) = ^^0(1 + '^i?^) = '-^1,7(7^(1 + '-^R.t) is the (time varying) Larmor 
frequency of the particle {cul is the usual cyclotron frequency). From the con- 
dition eq. (7) the i-th particle of the ensemble emits its synchrotron burst at 
the time (we retain here for obvious reasons only the positive solutions of eq. 

(7)) 

^. = ^[(l + 4^«.)^/^-l]- (8) 



From eq. (8) it is easy to derive the Lorentz factor of the particle at the 
moment of the burst: 



Eq. (9) will be used to rewrite the synchrotron spectrum as the spectrum of 
an ensemble of particles each characterized by a Lorentz factor 7j and a phase 
ai. 

In paper I we have discussed the differences in synchrotron emission produced 
by an ensemble of particles in a coherent and an incoherent configuration 
of phases. In the present paper we will only be concerned with the case of 
incoherent radiation, so that the energy radiated per unit frequency and unit 
solid angle in the directions perpendicular and parallel to the magnetic field 
can be written respectively as follows: 
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and 



dujd^l Aire 




(11) 



where we defined the usual quantities [3,4]: 



Ok = {i + eV,y/' 



,2 ki 



(12) 



and 9 = Ti /2 — 9z with 9^ the zenith angle of the versor n that points to the 
observer. 

The quantity that is most interesting from a physical point of view is the power 
(per unit frequency) radiated by the particles. To obtain this quantity, as we 
have discussed in paper I, one has to divide the energy per unit frequency over 
the observation time Tobs- The total power radiated per unit frequency by the 
Z particles in a period Tohs i^U 



where Xk = uj/{^ujliI)- 

Consider now the summation over the phases a^, we may consider these phases 
homogeneously distributed between (0, om), where au is the maximum value 
of a such that particles can illuminate the observer during the observation 
time Tobs [recall eq. (8)]: 



Introducing the phase density p{a) = Z/{27i), we may pass from a discrete 
sum to an integral by the substitution 



In conclusion, the power radiated per unit frequency by an ensemble of Z 
particles (all with the same initial Lorentz factor) in the backreaction regime 

^ Here we have used the standard computation of the solid angle integration. 




(13) 



= t^L7o ^Tobs{Tobs^^R + 2) . 



(14) 




(15) 
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will be 



dP V3q^ Z 1 f u! /■ 

a; (a) 

The spectrum found in eq. (16) has some general features that it is worth 
to investigate. First, let us check that in the limit Tobs <^ c^r^ the spectrum 
in eq. (16) converges to the standard synchrotron spectrum. This condition 
guarantees that for small observation times there is no appreciable difference 
between the well known synchrotron spectrum and the predicted one. It is 
useful to introduce the variable A = AujRTobsii^RTobs + 2) and rewrite the 
integral over a as an integral over y — 4^q;, so that the power radiated per 
unit frequency reads 



dP ^/3q^ Z (jjo 



du 8c 27r Jl + \K 



- / dyF{y), (17) 

1 n 



where we defined 



F{y) = livXv) I diK„^{i) . 

Now we compute the limit A of eq. (17) [note that this corresponds to 
evaluate the limit for cojiTobs ^ 1] : 

A 

lim , / dyF(y) = — F(0). (18) 



Substituting this expression in (16) the standard synchrotron power per unit 
frequency by an ensemble of Z particles all with the same Lorentz factor is 
readily recovered: 

Xo 



where Xq — u!/{3u!l^q). 

It is now useful, mainly for practical purposes, to derive analytical approx- 
imations for the spectrum of the radiation. We assume first to be in the 
regime where there are severe modifications due to the backreaction, that 
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is {ujRTobs > 1)- To perform the calculation we will adopt the following rough 
approximation of the Bessel function: 



X 




< x < 1 

x> I 



Although certainly not sophisticated, this approximation allows to treat an- 
alytically expressions that would otherwise be only of numerical access. The 
condition x{q.) < 1, using equation (9) and the definition of x, implies that 



a < 70 — ^ — 



7o - 



1/2 



(20) 



Therefore, to use the approximation of the Bessel function, we have to perform 
the integration over a between (0, a) where a — mm{aM, C(o}- The condition 
ao = ckm determines the frequency where there is a change in the slope of 
the spectrum of the emitted radiation. This identifies two frequency regimes, 
the low frequency one (a = Om) and the high frequency one (a — ao)- The 
separation between the two regimes occurs at 

u = 3a;i7oV(r„6.) (21) 



where 



Yin,,) 



'i^RTobsi^RTohs + 2) 

and, using the condition of strong backreaction ujrToi,s > 1, we may approxi 
mate 

1 



{ujRTohsf 



(22) 



In the high frequency regime ?)UJ l'~^IY (Tobs) < < SlulIq the integration over 
a has to be performed in the range (0, ao)- 

Within the approximations adopted here, we get the following expression for 
the radiated power: 

• Low frequency 
u < 3uL-fiY{Tobs) 
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High frequency 



dw c 27:5 \3j (jRTobs \^oul1o J 

Eqs. (23,24) are very interesting and allow a simple interpretation of the back- 
reaction regime. In the low frequency regime the spectrum of the synchrotron 
radiation has the same slope as in the standard case (uj^^^); on the other hand, 
moving to high frequency, the spectrum has a new power law behaviour of the 
type a;~^/^. Moreover, in the backrcaction regime, the maximum of the spec- 
trum is located at a frequency that depends on the observation time according 
to the following expression: 

u^^in^s) = ^LlMTo,s) ^ T^^^jA-, = , ■ (25) 

['^Rj-obs) y'^K-l-obs) 



In other words the maximum gradually moves toward lower frequencies when 
the observation time is increased and its position does not depend on 70: the 
integrated power becomes gradually richer in its low frequency component. 

We plotted our results for the radiation spectra in fig. 2 for three 
indicated: 1) standard case; 2) A = 10^ and 3) A = 10^. Larger values of A 
correspond to higher levels of backreaction, so that the spectra are gradually 
peaked at lower frequencies when the observation time increases. For the two 
cases of backreaction, we also plot our analytical approximations for the low 
and high frequency regimes. At low frequency the agreement with the detailed 
calculation is excellent. At higher frequencies clearly the rough approximation 
adopted for the Bessel functions gives a poorer but still acceptable agreement, 
very useful for practical estimates. In particular the slopes predicted by the 
analytical calculations reproduce very well the ones obtained in the detailed 
calculations. 

Note that the plot in fig. 2 is made in such a way that can be applied equally 
well to electrons and protons as radiating particles and for any value of the 
magnetic field. All these parameters in fact enter the definition of a;^,. 



4 Synchrotron radiation from particles with a power law spectrum 

We assume here to have a spectrum of charged particles in the form 

A^(7o) = A^o7o"" (26) 
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Fig. 2. Power spectrum of the radiation emitted by monoenergetic particles for the 
standard case ( no backreaction ) and for the cases A = 10^ and A = 10^ (The dashed 
lines represent the results of the analytical approximation illustrated in the text). 

where p > 1 is a spectral index. For most astrophysical applications this is 
the relevant case, therefore we explore here in detail what are our predictions 
for the spectrum of the synchrotron emission. Note that A^(7o) is here what is 
usually called the equilibrium spectrum of the radiating particles, as derived 
from a transport equation, usually including energy losses. This spectrum 
is typically steeper than the injection spectrum. The equilibrium spectrum 
is here taken to be time independent, meaning that there is a continuous 
replenishment of the particles at all energies, requiring the source to be active 
for the all duration of the phenomenon. The case of sources bursting on time 
scales much shorter than the observation time will not be considered here, but 
can be easily recovered by generalizing the discussion below. 

Starting from eq. (16) we can generalize the calculations of the previous section 
simply by the substitution Z — > iV(7o) and introducing an integration over 
the energy spectrum of the particles. The power radiated will be 

where 7^,^ and 

'Imax s-^s the minimum and maximum values allowed for the 
Lorentz factor, fixed by the particular physical system under consideration. 
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Using the same approximation of the Bessel function introdTiccd in the pre- 
vious section we will work out from equation (27) some interesting analytical 
results. Let us first assume that, fixing Tq^^, for any 70 inside the interval 
{lmin,1max) the system is in the backreaction regime uj^Tobs > 1- Using the 
high and low frequency results of the previous section (cfr. eqs. (24) and (23) 
we may calculate the integral over 70. 

Let us start from the high frequency regime, ?iU}l'~^IY {Tabs) < u) < 3a;x,7Q. In 
this case, fixing the frequency uj from the condition uj < 'iujLlo obtains 

therefore, for consistency, the integral over 70 has to be calculated between 7 
and 7maa:- this case (recalling eq. (24)) the power radiated per unit frequency 
is 

^^^^ 
7 



where A is a numerical factor defined as 

iVn 



27rc 5 



Assuming that '^max ^ 7 we may consider the integral in (28) from (7, 00) 
obtaining the high frequency power emitted per unit frequency 



J high ^KUbs P - 1 



Let us now evaluate the power emitted in the low frequency regime uj < 
2>u}Ly{Tobs)- Recalling that we are considering the case in which the backre- 
action is always important by construction, for any 70 G (cim.in,lmax)-i one can 
approximate Y{Tobs) — ^/ {yjRTobsf (cfr. eq. (22)) and therefore 

{uJKTobsY ' 

In this case there is no condition on the lower extreme in the integral over 
7o, so that, assuming ^min ^ Imaxi the low frequency power emitted per unit 
frequency can be written as 

(g_^^^...(..r^)-^(3^) . (30) 
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It is instructive to study the frequency corresponding to the change from the 
low to the high frequency regime. As we have pointed out before [eq. (22)], 
this frequency depends on the observation time according to the following 
expression: 

uj^ = ?,ujLllY{To,s) ^ ■ (31) 



This expression should be compared with the critical frequency defined within 
the standard theory of synchrotron emission for the particles with the mini- 
mum Lorentz factor, ujsyn — ^^LJmin- Since we assumed that, for any 70, we 
are in the backreaction regime {coK^rninTobs > 1), we have 

therefore, the change in the behaviour of the power radiated per unit frequency 
occurs in the low frequency regime. 

We are now ready to discuss a more realistic situation, for which there is a 
value 7c with ^ymin < 7c < Imax that divides, for any fixed observation time, 
the standard regime from the backreaction one. This value of % is simply 
given by 

_ 1 

In the situation for which 7c G {"^mmi Imax) we can divide the spectrum in two 
regions: one for frequencies uo < SuJi'-f'^, in which backreaction is not important, 
and the other, for u > Sul^"^, in which the backreaction effects need to be 
included. 

In the situation in which 7c £ {■Jmin, Jmax) the spectrum radiated by particles 
is 



At high frequencies, cu > 80;^ 7^: 

dP 




back 

(32) 



At low frequencies, uj < SujlJ, 



dP fdPV' (dPV"'^ 




duj 

where the first term is the standard synchrotron spectrum. 



(33) 
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The standard synchrotron spectrum (integrated between ^min and 7c) has the 
following form: 



do; y ° c 27r 



37 3p - 1 



p-1 

2 



3a;z, 



7c' 



-p+i 



(34) 



while the backreaction dominated spectrum is described by eqs. (29) and (30) 
with 7^j„ = 7c. Let us now compare the two components of the low frequency 
spectrum. For this purpose we evaluate explicitly eq. (33) 



C ■ 



3p — 1 y^ujL 




hack 



low 




£-1 

2 



obs) 



3a; r 



3 1 



5 p — 1 3p — 1 



(35) 



where C is a numerical factor 

c 271 V3 
and we have used the relation 7c = 1/ {utxTobs) ■ 



It is easy to see that at frequencies cu <S cUcr the low frequency spectrum is 
well described by the standard synchrotron emission (i.e. the second term in 
(35) is always negligible). Moreover, this is also rigorously true in the special 
case p = 2, as it is easy to show by explicitly calculating the second term in 
eq. (35). 

In conclusion, in the most interesting case in which radiation reaction is effi- 
cient only for 7 > 7c > 7mm) the power radiated by the ensemble of particles 
with a fixed (power law) distribution of Lorentz factors is described at low fre- 
quency {cu < Sujl^'^) by the standard synchrotron spectrum (34) and at high 
frequency {cu > Scul^'^) by the backreaction modified spectrum in eq. (29). 

In fig. 3 we plot for illustrative purposes the results of a specific case: high 
energy protons with a spectrum oc E~'^ and high energy cutoff at 10^^ eV 
radiating by synchrotron emission in magnetic fields of order 100 Gauss (solid 
lines) and 1000 Gauss (dashed lines). The curves are obtained for three choices 
of the observation time Tots- The rightmost curves in both cases correspond 
to the standard synchrotron emission (small observation time) while moving 
leftward we plot the cases To^^ = 10^ s and Tg^s — 10^ s. The effect of the 
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Fig. 3. Power spectrum of the radiation emitted by protons with a power law spec- 
trum, N{^q) = Nq^q^ for the standard case (T^bs ~^ ^) o-iT'd for T^hg = 10^ s and 
Tobs = 10*^ s (moving leftward). The two sets of curves are for B = 100 Gauss (solid 
lines) and B = 1000 Gauss (dashed lines). 

increasing backrcaction is clear from these curves. This specific case has been 
chosen because the situation is similar to the one proposed in [5,6] to explain 
the TeV emission from BL Lac as synchrotron emission of ultra high energy 
protons (in these papers the standard synchrotron emission was adopted). 



5 Conclusions 



We calculated the spectrum of the synchrotron emission from a system of 
charged particles in the most general case in which the observation time is 
arbitrary compared to the time for the energy losses of the particles. We first 
calculated the effect for the case of a particle with Lorentz factor 79. 

Our findings on the monoenergetic case can be summarized as follows: 

i) There is a critical frequency Ucr such that the spectrum of the radiation 
at a; < Ucr is the usual synchrotron spectrum oc uj^^^, while for uj > Ucr the 
backreaction affects the spectrum changing it to oc o;"^/^. 

ii) The critical frequency ujcr depends on the observation time but it turns 
out to be independent of the Lorentz factor of the radiating particles 70. The 
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dependence on these parameters is as found in eq. (25), so that increasing the 
observation time, power is moved to gradually lower frequencies. 

in) The standard limit is recovered when the observation time is much smaller 
than the time for losses unTobs ^ 1- 

The more realistic case investigated in this paper is that of a power law spec- 
trum of particles NIjq) oc 7^^^. In general what happens is that there may 
be a critical Lorentz factor 7c, such that particles with 70 > 7c are affected 
by the backreaction while the particles with 70 < 7c behave in the standard 
way. The spectrum of the radiation in this case is a superposition of different 
components. Our findings can be summarized as follows: 

a) the spectrum of the radiation can be divided into a low frequency one 
and a high frequency one, with the separation occurring at the frequency 

OJcr — "i^L/i^KTohsf ■ 

b) the particles with Lorentz factor 70 < 7c radiate the standard synchrotron 
radiation whose spectrum is oc u'^'^^^^. 

c) the particles with 70 > 7c do radiate in regime of backreaction and affect 
both the low frequency and the high frequency regime. 

d) the low frequency spectrum radiated by particles with 70 > 7c is oc u^^^ 
for uj <^ ujcr, but it is always dominated, in the same frequency range, by the 
standard synchrotron radiation. 

e) the high frequency radiation radiated by particles with 70 > 7c is oc ijj^'^I'^. 
Therefore it represents a suppression of the radiation compared to the case 
of standard synchrotron radiation. Note that increasing the observation time, 
while the slopes at low and high frequency remain unchanged, the boundary 
between the two regimes moves toward lower frequencies, so that, as a conse- 
quence, the height of the spectrum at high frequencies becomes increasingly 
lower. 

f) At fixed observation time (which is obviously decided by the observer) 
the critical frequency only depends on the magnetic field in the production 
region. This is a very important point: in the standard synchrotron emission, 
it is in general not possible to extract the magnetic field from a measurement 
of the synchrotron flux because there is degeneracy between the number of 
radiating particles and the strength of the magnetic field. In the backreaction 
case however, the position of the change in slope uniquely defines the magnetic 
field, so that the measurement can in principle be used to directly infer the 
strength of the magnetic field. 

The crucial question is whether there are situations in which the conditions 
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for the backreaction to be relevant arc fulfilled. The answer can be found in 
fig. 1, where we plotted the regions of interest for both electrons and protons 
in the plane B — E. At each magnetic field there corresponds a range of 
energies for which the backreaction is important. We immediately see that 
for conditions typical in the Galaxy, B ~ IjiG only electrons with E ^ 10^ 
GeV feel the effects of the backreaction. On the other hand in magnetic fields 
which arc typical of neutron stars, B ~ 10^° — lO"*^^ Gauss, electrons with 
energies in excess of MeV-GeV already need to be accounted for in the frame 
of a backreaction approach. In the same environment, protons with energies 
larger than 10^ — 10^ GeV also radiate in a regime in which backreaction is 
important. 

Some applications of the calculations and results illustrated in this paper will 
be presented in forthcoming papers. They include: the synchrotron emission 
from ultra-high energy electron-positron pairs generated as a result of the 
decay of super-heavy relics [7] in the Galaxy [8] or in the decay of the Z° 
resonance produced in uu annihilation (Z-burst model [9,10]). In this case we 
have magnetic fields B ~ l/xG and energies in excess of 10^° GeV, so that the 
synchrotron radiation is strongly affected by the backreaction. 

A more conventional application concerns the TeV gamma ray emission from 
BL Lac objects. In [5,6] a proton synchrotron model was proposed in which 
the TeV emission is the result of the synchrotron emission of ultra high energy 
protons {E ~ 10^^ eV) in a magnetic field of order ~ 100 Gauss. We estimate 
that with these parameters and for observation times of the order of the du- 
ration of the observed flares (a few hours) the backreaction affects visibly the 
TeV gamma ray spectra [11]. 
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